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Graphene samples can have a very high carrier mobility if influences from the sub¬ 
strate and the environment are minimized. Embedding a graphene sheet into a het¬ 
erostructure with hexagonal boron nitride (hBN) on both sides was shown to be a par¬ 
ticularly efficient way of achieving a high bulk mobility [!]• Nanopatterning graphene 
can add extra damage and drastically reduce sample mobility by edge disorder nni. 
Preparing etched graphene nanostructures on top of an hBN substrate instead of Si02 
is no remedy, as transport characteristics are still dominated by edge roughness [^. 
Here we show that etching fully encapsulated graphene on the nanoscale is more gentle 
and the high mobility can be preserved. To this end, we prepared graphene antidot 
lattices [6] where we observe magnetotransport features stemming from ballistic trans¬ 
port. Due to the short lattice period in our samples we can also explore the boundary 
between the classical and the quantum transport regime. 


In single layer graphene the charge carriers are com¬ 
pletely exposed to the environment, which limits their 
mobility. Placing graphene on hexagonal boron nitride 
(hBN) was shown to improve the carrier mobility [7], al¬ 
lowing the observation of ballistic transport or the frac¬ 
tional quantum Hall effect in bulk graphene [5] . Recently, 
a dry stacking technique was introduced, which allows 
complete encapsulation of graphene into layers of hBN 
and excludes any contamination from process chemicals 
such as electron beam resist [T]. To obtain graphene 
nanodevices, chemically prepared graphene nanostruc¬ 
tures [MI] are a potential route for certain applications, 
however, the high flexibility of a top down patterning 
approach is extremely desirable. Graphene antidot lat¬ 
tices can help circumventing the problem of the miss¬ 
ing band gap in transistor applications m, and were 
even predicted to serve as the technological basis for spin 
qubits [13]. Clearly, for advanced graphene nanodevices, 
not only the bulk mobility has to be improved, but the 
nanopatterning has to be optimised. 

Here we present experiments on graphene antidot lat¬ 
tices [3 m ng etched into hBN/graphene/hBN het¬ 
erostructures with lattice periods going down to a = 50 
nm. Magnetotransport on those samples shows commen- 
surability features stemming from ballistic orbits around 
one or several antidots. This allows us to prove that the 
high carrier mobility is preserved in the nanopatterning 
step even though the zero field resistance is dominated by 
scattering on the artihcial nanopattern, giving an appar¬ 
ent reduction of the mobility. The small feature size of 
our samples also allows us to approach the region where 
the classical picture of cyclotron orbits no longer applies. 
This classical to quantum crossover is governed by the 
ratio between the Fermi wavelength of the carriers 
and the dimensions of the nanopattern. 

To obtain embedded graphene samples. 


hBN/graphene/hBN stacks were prepared using the 
dry stacking technique, patterned into Hall bar shape, 
and contacted using Cr/Au [T]. In hBN/graphene/hBN 
samples prepared by this method, we routinely obtained 
carrier mobilities in excess of /r = 100 000 cm^/Vs, 
showing all integer quantum Hall states starting from a 
few Tesla. In one sample without antidots and a mobility 
of /r = 300 000 cm^/Vs we also observed the fractional 
quantum Hall effect at T = 1.4 K. This shows that our 
fabrication procedure is mature and consistently yields 
high sample qualities. The samples presented in this 
study did not show any signs of a moire superlattice 
[mu?]. Afterwards, an antidot lattice was patterned. 
(For fabrication details, see Methods section). Fig. 
1 shows an optical micrograph of a finished sample 
and a scanning electron micrograph of a sample after 
measuring. A sketch of the antidot lattice, etched into 
the stack is also shown. The antidot lattice period a 
was varied between 50 nm and 250 nm. The antidot 
diameter d was lithographically defined to be about 40 
nm, but due to the conical etching profile, the actual 
diameter in the graphene plane is smaller. Using SEM 
inspection, we estimate it to be about 25... 30 nm. 

In Fig. 2, we show data for a sample with a lattice 
period of a = 200 nm. From the gate response of the con¬ 
ductivity at a magnetic field B = 0, shown in Fig. 2a, we 
calculate an apparent field effect mobility of p, = 35 000 
cm^/Vs. At a carrier density ns = 2.3 x 10^^ cm“^ 
this corresponds to an apparent mean free path of about 
^mfp = = 620 nm. We estimate the intrinsic 

mean free path to be about 1400 nm [TS] (see Supplemen¬ 
tary Information). Magnetotransport traces of this de¬ 
vice (see Fig. 2b) show pronounced peaks at held values 
where the cyclotron diameter 2Rc = is com¬ 

mensurate to the square antidot lattice. The peak be¬ 
longing to 2Rc = a, the fundamental antidot peak, is 
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most pronounced. Additional peaks appearing at lower 
fields correspond to orbits encircling 2, and 4 antidots [6] 
(see Fig. Id), confirming a mean free path which spans 
several lattice periods. While in a simple picture only the 
unperturbed orbits encircling the antidots are responsi¬ 
ble for the magnetotransport features, a more detailed 
analysis based on the Kubo formula shows that veloc¬ 
ity correlations in the chaotic trajectories, which occupy 
the largest part of the phase space, result in the magne¬ 
toresistance peaks [niiin]. Most of the orbits therefore 
hit the antidot edges several times within a mean free 
path. Hence, the visibility of the antidot peaks not only 
proves a high bulk mobility, but also shows that scatter¬ 
ing at the edges does not cut off the trajectories and we 
can conclude that the high carrier mobility also survives 
after nanopatterning. 

At higher fields, the cyclotron diameter 2Rc is reduced 
below the neck width a — d in between the antidots. 
We can observe Shubnikov-de Haas oscillations, eventu¬ 
ally resulting in a well-defined quantum Hall effect. At 
H = 14 T we clearly observe the ly = 1 plateau, which 
again shows the high sample quality (Fig. 2c). We eval¬ 
uated the carrier density dependence of the magnetore¬ 
sistance peaks corresponding to orbits around 1, 2 and 
4 (Fig. 2d) and found that the peaks were always well 
described by a square root dependence of the cyclotron 
diameter on the carrier density down to ns = 3.2 x 10^^ 
cm“^. Quantitatively, we confirmed the formula for the 
cyclotron diameter for graphene given above, which con¬ 
tains spin and valley degeneracy. 

Fig. 3a shows the magnetoresistance of a sample with 
a = 100 nm at ns = 2.8 x 10^^ cm“^. The appar¬ 
ent Hall mobility at this density is about n = 8 000 
cm^/Vs. Again, scattering at the antidot potential lim¬ 
its the apparent mobility |18j . but the intrinsic mobility 
is higher as we clearly observe magnetoresistance peaks 
for n = 1, 2,4 antidots, and a fourth peak at lower fields 
is weakly visible. Ishizaka and Ando studied how the 
visibility of the higher order antidot peaks depends on 
the mobility [2^. From their data, we estimate that the 
intrinsic mean free path must be at least 400 nm, well 
in excess of the apparent mean free path of 160 nm (see 
Supplementary Information). 

The good visibility of the n = 2 peak confirms the 
small aspect ratio dja < 0.3 [5D], in agreement with our 
SEM analysis and also with the onset of the well-defined 
Shubnikov-de Haas oscillations in our magnetotransport 
data. All these approaches give an antidot diameter of 
d = 25 ... 30 nm. 

In experiments in GaAs based antidot lattices it was 
found that due to depletion at the antidot boundaries, 
the potential can be very soft and small lattice periods 
are hard to realize. In our case the data compares well 
to hard-wall potential lattices in GaAs, which could be 
realized in GaAs only at much larger lattice periods [B]. 
We also compared data for similar carrier densities in 


the electron and hole regime in Fig. 3b and found the 
graphs to be virtually identical. This proves that there is 
no edge doping at the antidot boundaries, which would 
have led to different potential shapes in the electron and 
hole regime due to Fermi level pinning at the edges. 

Now let us discuss the transition between the quantum 
and the classical transport regime. In GaAs-based het¬ 
erostructures, the smallest lattice period realised so far 
was o = 80 nm, and required critical tuning of the etch 
depth [^. In contrast, due to the lack of a depletion 
region in graphene the fabrication of samples with a very 
small lattice period is less critical, and the carrier den¬ 
sity is widely tunable. Also, due to valley degeneracy, 
the Fermi wavelength in graphene, Xp = is a fac¬ 

tor of y/2 larger than in GaAs based 2DEGs at the same 
carrier density. Thus we can explore the transition from 
the semi-classical to the quantum regime [22], where a 
description in terms of classical orbits is no longer jus¬ 
tified. In the samples with a < 100 nm we are able 
to study this transition. Fig. 4a shows the disappear¬ 
ance of the main antidot peak in a sample with a = 75 
nm as the carrier density is lowered, making Xp longer. 
We find that this peak is only visible at densities above 
ns = 4.3 X 10^^ cm“^, corresponding to Xp = 54 nm. 
Also, in two samples with a = 100 nm, we observe that 
the main antidot peak becomes visible for densities larger 
than 715 = 2.2x10^^ cm“^, which corresponds to Xp = 75 
nm. In a sample with a = 50 nm, we observed a weak 
antidot peak only at ns = 2.5 x 10^^ cm“^(A_F = 22 
nm). To be in the classical limit of a quantum system, 
the Fermi wavelength must satisfy a condition ^ I 
[23] ■ where I is a typical dimension of the system. In our 
case, the neck width a — d of the constriction between the 
antidots is the shortest length scale in the problem, and 
we find that when Xp ~ a — d the classical regime sets in 
and the antidot peak becomes visible. 

The fact that the antidot peaks disappear at low den¬ 
sities can be either due to a limited mean free path or 
the breakdown of the classical picture. In Fig. 2d (lat¬ 
tice period a = 200 nm) all the antidot peaks disappear 
at roughly the same magnetic field, B « 0.5 T (where 
/iH exceeds some constant), but different carrier density. 
This behaviour is clearly governed by a limited mean free 
path. In contrast, in the sample of Fig. 4a (a = 100 nm), 
we find that the classical features at both H « 1 T and 
H « 2.5 T disappear at the same carrier densities, mak¬ 
ing a Ai?-driven scenario more realistic. 

Finally, at low densities, we can observe a weak local¬ 
ization (WL) feature at low temperatures: a peak in the 
magnetoresistance at B = 0 (see Fig. 4b). Using a stan¬ 
dard analysis for WL in graphene |24| that we employed 
in earlier work on graphene antidot lattices on Si 02 
we extracted the phase coherence length L^. For the 
sample with a = 100 nm (same as in Fig. 3a) we found 
it to be between 120 nm and 300 nm (see Fig. 4c). It 
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clearly exceeds the lattice period, unlike in graphene an¬ 
tidot samples on Si 02 where was significantly below a 
[14] . We therefore again conclude that nanopatterning of 
embedded graphene leads to greatly reduced scattering 
at the sample edges. 

In summary, we prepared antidot lattices in stacks of 
hBN/graphene/hBN and observed well-developed com- 
mensurability features in samples with lattice periods 
from a = 50 nm to a = 250 nm. This shows that the 
etching procedure preserves the high sample quality. In 
the short-period graphene samples, we could observe the 
disappearance of classical features when the Fermi wave¬ 
length A F exceeds a — d, marking a classical to quantum 
transition. Our experiments therefore pave the way for 
well-controlled graphene based nanodevices. 

METHODS 

Single crystalline hexagonal boron nitride fhBNipSl 
im was exfoliated onto a stack of PMGI and PMMA 
polymers spin coated on an oxidised Si-Wafer [7]. Suit¬ 
able hBN flakes, serving later as the top hBN layer, were 
located in an optical microscope by using different band¬ 
pass hlters. The PMGI sacrificial layer was dissolved in 
photoresist developer and DI water, leaving the PMMA 
with the hBN floating on DI water. Then the PMMA 
film was transferred to a microscope slide into which a 
hole had been cut [7]. Single layer graphene was exfoli¬ 
ated from HOPG (Momentive Performance, ZYA grade) 
onto oxidised Si, and picked up by van-der-Waals inter¬ 
action using the first hBN flake [1]. Using a home-made 
setup in an optical microscope, this stack was transferred 
to a second hBN flake, residing on an oxidised Si sub¬ 
strate with prepatterned markers. Subsequently, the fin¬ 
ished stacks were annealed in forming gas flow at 320° G 
for several hours. The hBN/graphene/hBN heterostruc¬ 
ture was then patterned into Hall bar shape using elec¬ 
tron beam lithography (EBL) and CHF 3 /O 2 (40 sccm /6 
seem, 60 Watt power) based reactive ion etching (RIE) 
[T]. Cr/Au side contacts were defined with EBL and 
deposited by thermal evaporation and lift-off after brief 
oxygen plasma cleaning of the contact areas [1]. Finally, 
the antidot lattice was defined in a separate EBL and 
RIE step. Samples were glued into chip carriers with sil¬ 
ver filled epoxy to contact the back gate, wire-bonded 
and measured in a helium cryostat with variable temper¬ 
ature insert, using low frequency lock in techniques with 
a bias current of 10 nA. 
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FIGURES 


FIG. 1. a, Optical micrograph of a finished graphene Hall bar 
etched out of an hBN/graphene/hBN heterostructure, and 
contacted with Cr/Au leads. Scale bar length 5 fim. b, False- 
color scanning electron micrograph of a sample with lattice 
period a = 100 nm. Here, the heterostructure is shaded in 
green, the Cr/Au contacts yellow, and the Si/Si02 substrate 
violet. Scale bar length 500 nm. c. Sketch of the antidot lat¬ 
tice in an hBN/graphene/hBN heterostructure. The antidot 
lattice period a ranges from 50 to 250 nm, and the antidot 
diameter d is about 25... 30 nm. d. The most prominent 
cyclotron orbits fitting into the lattice, giving rise to magne¬ 
toresistance peaks. 
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FIG. 2. a, Gate dependence of the sheet conductivity of a 
sample with a = 200 nm. The linear fit (red line) gives an ap¬ 
parent mobility of /i = 35 000 cm^/Vs. b, Magnetoresistance 
(black) and Hall resistance (red). The arrows correspond to 
the expected magnetic field positions of the orbits sketched 
in Fig Id. The fine structure in is not noise, but phase- 
coherent oscillations that disappear at higher temper¬ 

atures. c, Gate dependence of Rxx and Rxy at B = 14 T, 
showing a clear u = 1 quantum Hall plateau. This feature 
is only observed in high-mobility graphene devices, d, The 
magnetic field positions of the three antidot peaks scale with 
the square root of the carrier density, confirming the classical 
nature of those peaks. 
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FIG. 3. a, Magnetoresistance and Hall resistance data taken 
on a sample with a = 100 nm. The three well-defined antidot 
peaks correspond to orbits around 1, 2 and 4 antidots, b, Rxx 
taken at similar electron and hole density. Here, T = 80 K to 
show the classical features more clearly. There is virtually no 
difference between those graphs, proving that the potential 
profile is the same for electrons and holes. 
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FIG. 4. a, Rxx data of a sample with a = 75 nm taken at 
very low densities, at the transition into the regime of clas¬ 
sical transport. The densities ns are given in units of 10^^ 
cm“^, shown next to the corresponding graphs. The expected 
position of the main antidot peak is marked with a triangle 
for each density. A higher order antidot peak is also visi¬ 
ble at lower magnetic fields. As the carrier density is low¬ 
ered, the antidot peaks disappear. Inset: Sketch of the Fermi 
wavelength corresponding to the red graph, b. Weak local¬ 
ization (WL) peak in the sample with a = 100 nm, taken at 
ns = 1.3 X 10^^ cm“^. The antidot peak is not visible at this 
low density, the big peaks at B = ±1.4 T are a Shubnikov-de 
Haas oscillation, c. Phase coherence length taken from the 
weak localization fits of a sample with a = 100 nm at various 
low densities and T = 1.4 K. The phase coherence length ex¬ 
ceeds the lattice period, showing that the etched boundaries 
do not lead to severe phase-breaking. 









